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We consider the stationary incompressible Navier-Stokes equation in the half-plane 
with inhomogeneous boundary condition. We prove existence of strong solutions for 
boundary data close to any Jeffery-Hamel solution with small flux evaluated on the 
boundary. The perturbation of the Jeffery-Hamel solution on the boundary has to 
satisfy a nonlinear compatibility condition which corresponds to the integral of the 
velocity field on the boundary. The first component of this integral is the flux which 
is an invariant quantity, but the second, called the asymmetry, is not invariant, which 
leads to one compatibility condition. Finally, we prove existence of weak solutions, 
as well as weak-strong uniqueness for small data. 
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1. Introduction 

The stationary and incompressible Navier-Stokes equations in the half-plane 

^ = {{x,y) e y > l} 


are 


Au —'Vp = u •'Vu , V-u = 0, 
'*^\dn = lim u = 0 , 

|aj|^oo 


( 1 ) 


where u* is a boundary condition. Due to the incompressibility of the fluid, the flux is an 
invariant quantity, 

u*-n= v{x,y)dx, 

Jdn Jr 

for all y > 1, where u = {u,v) and n = (0,1) is the normal vector to the half-plane. This 
problem (see figure la) presents three difficulties: D is a two-dimensional unbounded domain, 
the boundary of 5D is unbounded, and the boundary data are not zero. There is not much 
previous work on this problem, but some authors have treated related problems. Concerning the 
half-plane problem, Heywood (1976, §5) proves the uniqueness of solutions for the steady Stokes 
equation and the time-dependent Navier-Stokes equation. The so called Leray’s problem, which 
consists of a finite number of outlets connected to a compact domain, has been studied in detail 
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Figure 1: (a) The domain ri we consider is the half plane dehned by x £ M and y > 1; (b) an 
aperture domain; (c) a channel connected to a half-plane. 


by Amick (1977, 1978, 1980) and several other authors, but the resolvability for large fluxes 
is still an open problem. Fraenkel (1962, 1963) provides a formal asymptotic expansion of the 
stream function in case of a curved channel by starting with the Jeffery-Hamel solution (Jeffery, 
1915; Hamel, 1917) for the hrst order. The case of paraboloidal outlets was first treated by 
Nazarov &: Pileckas (1998), and then more recently by Kaulakyte &: Pileckas (2012); Kaulakyte 
(2013). Another important class of similar problems are the aperture domains, introduced by 
Heywood (1976), as shown in figure lb. The linear approximation was studied in any dimension 
by Farwig (1996); Farwig &: Sohr (1996). The three-dimensional case was treated by Borchers 
&: Pileckas (1992), as well as other authors. For the two-dimensional nonlinear problem, Galdi 
et al. (1995) proved that the velocity tends to zero in the L^-norm for arbitrary values of the 
flux. For small fluxes, Galdi et al. (1996) and Nazarov (1996) show that the asymptotic behavior 
is given by a Jeffery-Hamel solution but only if the problem is symmetry with respect to the 
y-axis. The asymptotic behavior of the two-dimensional aperture problem in the nonsymmetric 
case is still open. Finally, Nazarov et al. (2001, 2002) considered a straight channel connected 
to a half-plane (see hgure Ic), and looked under which conditions the asymptotic behavior is 
given by a Jeffery-Hamel flow in the half-plane and by the Poiseuille flow in the channel. Theses 
conditions are described in details later on. On a more applied side, the bifurcation properties 
and the stability of the Jeffery-Hamel flows have retained the attention of many authors (Moffatt 
Duffy, 1980; Sobey Sz Drazin, 1986; Banks et al., 1988; Uribe et al., 1997; Drazin, 1999). 
Jeffery-Hamel flows play an important role in the asymptotic behavior of flows carrying flux. 
They own their name to the work of Jeffery (1915); Hamel (1917), and are radial scale invariant 
solutions of the two-dimensional stationary incompressible Navier-Stokes equations 

Au — {u • V) u — Vp = 0 , V • It = 0 , 

in domains 

D = |(r sin(0), r cos(0)) £ : r > 0 and 9 £ (—/3, /3)| 

with 13 £ (0, |], satisfying the boundary condition 

'“laD\{o} ^ ^ • 

Explicitly, a Jeffery-Hamel solution itje is of the form 

itjH(r,6') = -f{e)er, 
r 

with / a solution of the nonlinear second order ordinary differential equation 

/" + /'+4/ = 2C, 
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with C G M, satisfying the boundary condition f{±(3) = 0. The constant C is related to the flux 

of the flow, 

r+P 

/ f{e)de. 

j-fi 

The Jeffery-Hamel solutions have been intensively studied (Rosenhead, 1940; Fraenkel, 1962; 
Banks et ai, 1988; Tutty, 1996), and, because some mathematical questions still remain open, 
there has been a regain of interest in recent years (Rivkind &: Solonnikov, 2000; Kerswell et al, 
2004; Putkaradze k, Vorobieff, 2006; Corless k Assefa, 2007). 

In what follows we are interested in the half-plane case, so we consider /3 = i.e., when the 

domain is the upper half plane = M x (0, oo). In Cartesian coordinates, the two components 
of the velocity of the Jeffery-Hamel solutions are 

Vm{x,y) =-fv , (2) 

y \yJ 

/(arctan s) 

/»<») = l +,2 ■ 

The Jeffery-Hamel solutions for /3 = f have a peculiar prop¬ 
erty: for small $ < 0 there is more than one solution. In 
fact, as shown in appendix §A, $ = 0 is a tri-critical bifur¬ 
cation point, see figure 2. For small $ > 0 the Jeffery-Hamel 
problem has a solution which is symmetric with respect 
to the y-axis. The solution also exists for small values of 
4> < 0, but when crossing from 4> > 0 to 4> < 0, an additional 
pair of asymmetric solutions (related to each other by 
a reflection with respect to the y-axis) appears. For <I> = 0, 
the Jeffery-Hamel solution is the zero function, and will be 
ignored in what follows. 

The central idea of the method we use to study (1), is to 
interpret the system as an evolution equation with y playing 
the role of time. The boundary data of the original problem 
then become the initial data for the resulting Cauchy problem. 

This allows discussing the “time” dependence of quantities 
like the flux 4) = J’jju(x,y)dx and the asymmetry A = J^u{x,y)dx in a natural setting. We 
assume for the moment sufficient decay for these integrals to make sense. As can be seen from 
(2), the flux and the asymmetry are invariants for a Jeffery-Hamel solution itj, i.e., they are 
independent of the time y, and therefore 



Figure 2: Existence of the Jeffery- 
Hamel flows for small values of 
the flux 4). For 4* > 0, there ex¬ 
ists one symmetric solutions, but 
for 4) < 0, also two additional 
asymmetric solutions exist. 


Um{x,y) = -fu f-) , 

y \yj 


with 


fu{s) = 


s /(arctan s) 


1 -I- s' 


2 ’ 


A= u%{x,l)dx= / u%{x,y)dx= lim / u%{x,y)dx, 

Jr Jr Jr 

= / v%{x,l)dx = / v%{x,y)dx = lim / v%{x,y)dx. 

Jr Jr y^oo 

The flux is an invariant of the Navier-Stokes equation (1), i.e., if u is a solution of the Navier- 
Stokes equation, then for y > 1, 


v{x, 1) dx 


v(x, y) dx = lim 

y^oo 


v{x,y) dx. 
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but the asymmetry is not an invariant, so typically, 


u{x, 1) dx / 


u(x, y) dx 7 ^ lim 

y—^oo 


u{x, y) dx. 


As we will see below, for solutions that are close but not equal to Jeffery-Hamel, the asymmetry is 
no more an invariant, and this fact is the main source of trouble for the construction of solutions. 

Jeffery-Hamel solutions are singular at the origin, and, in order to study this phenomenon, 
i.e. look for solutions which are close to Jeffery-Hamel flows, it is necessary to regularize 
the problem. For this purpose, given a Jeffery-Hamel solution uj in the upper half plane 
Z) = {(x, y) G : y > 0}, we restrict it to the domain H = {(x, y) G : y > 1}, and construct 
stationary solutions of the Navier-Stokes equations which are close to the Jeffery-Hamel flow by 
imposing boundary conditions of the form 




(3) 


with small and with zero flux. 


Uf, • n = 0 . 


I an 


Even when considering such boundary conditions, we are not able to perform a fixed point 
argument on the nonlinearity by inverting the Stokes problem. The main reason is that the flux 
is determined by the boundary condition, while the asymmetry is not. In order to adjust the 
asymmetry, we rewrite the boundary condition as 


u\ 


an 


— '“'SUn + (^’0) 


1 




e 2 ^ 


(4) 


where 


A = 


Ub , 


so that Ug has no asymmetry and no flux. 


U = 0. 


1 2 

The choice of e“ 2 ^ in (4) is for convenience later on, and we could have chosen instead any 
other smooth function of rapid decay. 

We will show the existence of strong solutions to the Navier-Stokes equation (1), with the 
boundary condition (4) for 4> small and Ug in a small ball by adjusting the parameter A. The 
main result is the following: 

Theorem 1. For all boundary condition of the form (4) with itj a Jeffery-Hamel solution with 
small enough flux and all Ug in a small enough neighborhood of zero in some function space, 
there exists a solution {u,p) of the Navier-Stokes equation in satisfying 


and 


lim y suplti — u%\ = 0 , 

Va:eK 7 


Vu G , yu G L°°{n ), 

u/y G L\n ), yVu G . 


(5) 


Moreover, if v is a weak solution (defined in (15)J of (1) and u a strong solution of (1) satisfying 
(5), then u = v, provided u* is small enough. 
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We now discuss more precisely the results of Nazarov et al. (2001, 2002), for the domain 
shown in figure Ic. We note that in this domain, the flux through the channel is not prescribed. 
They show that by requiring the asymptotic behavior to be an antisymmetric Jeffery-Hamel 
solution there exists a unique solution in some weighted space, and the flux is uniquely 
determined by the data. Conversely, by requiring that the asymptotic behavior is given by a 
symmetric Jeffery-Hamel solution the Navier-Stokes equation linearized around leads to 
a well-posed problem for < 0 and to an ill-posed one for $ > 0. So for $ < 0, there exists 
a unique solution for all small enough fluxes, but for $ > 0, the asymptotic behavior is still 
unknown. So we believe that in case <1> < 0, the Navier-Stokes equation in the half-plane (1) has 
a solution decaying like r~^ at inhnity whose asymptote is given by a Jeffery-Hamel solution, 
but in the case > 0, it is still not clear that the solution is in general bounded by r~^. 

The remainder of this paper is organized as follows. In section §2, we introduce the function 
spaces which we use for the mathematical formulation of the problem, and prove some basic 
bounds. In section §3, we rewrite the Stokes equation as a dynamical system, present the 
associated integral equations, and provide bounds on the solution of the Stokes system, so that 
in section §4, we can show the existence of strong solutions to the Navier-Stokes system. In 
section §5, we prove existence of weak solutions, and, finally, in section §6, we prove uniqueness 
of solutions for small data with a weak-strong uniqueness result. In the last part, we also present 
numerical simulations that show that the asymptotic behavior is most likely not given by the 
Jeffery-Hamel solution if > 0. In the appendix, we show the existence of symmetric and 
asymmetric Jeffery-Hamel solutions with small flux. 

2. Function spaces 

As explained in the introduction, our strategy of proof is to rewrite (1) as a dynamical system with 
y playing the role of time. This system is studied by taking the Fourier transform in the variable 
X, which transforms the system into a set of ordinary differential equations with respect to y. 
We now define the function spaces for the Fourier transforms of the velocity field, pressure field 
and the nonlinearity. The choice of spaces is motivated by the scaling property of the equations 
with respect to x and y when linearized around a Jeffery-Hamel solution. This setup turns out 
to be natural for the description of the asymptotic behavior of solutions close to Jeffery-Hamel 
flows. Similar function spaces were already used by Wittwer (2002), where the basic operations 
which are needed for the discussion of the Navier-Stokes equations were discussed. In particular, 
Wittwer (2002) shows basic bounds on the convolution with respect to the variable k, the Fourier 
conjugate variable of x, which is needed to implement the nonlinearities, and bounds on the 
convolution with the semigroup which is associated with the Stokes operator when viewed 

as a time evolution in y. Further properties and improved bounds have been proved by Hillairet 
Sz Wittwer (2009); Boeckle &: Wittwer (2012). 

Definition 2 (Fourier transform and convolution). For two functions / and g dehned almost 
everywhere in fl and which are in L^(M) for all y > 1, the inverse Fourier transform of / is 
defined by 

f{x,y) = T[f]{x,y) = [ e'^^/(A:, y) dfc , 

Jr 

and the convolution by 

{f* 9 ){k,y)= [ f{k- i,y)g{e,y)di. 

Jr 

We note that with these definitions, 

fg = J^[f*g] ■ 
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We now define two families of function spaces: the first one is for functions of k only which 
will be used for the boundary data, and the second one is for functions of k and y: 


Definition 3 (function spaces on tin). For a > 0 and q gM., let Aa,q be the Banach space of 
functions / G C{X,C) where 

9 > 0 , 

\{ 0 }, q<0, 


X = 


( 6 ) 


such that f{k) = f{—k) and such that the norm 

I/I 


/; Ac,q = sup 


X Vct,q 


with ya,q{k) = 


\k\^ 


1 + 


\a+q 


is finite. For a > 0 and (7 > 0, let Ta,q and Wa,q-i be the Banach space of functions in 
•^Q+i,min(o,<j-i) such that their respective norm 


Bl 

||/) "Toigll ^ ' 11 /) •^o+l,min(0,<7—1—i) || ) 

i=0 

is finite, where [gj denotes the integer part of q. 


Bl 

|/;w«,,|| = Ell^fc/;“^«+w-i-ilh 

i=0 


Remark 4. The parameter a captures the behavior of functions at infinity, which corresponds to 
the regularity in x in direct space. For example, if / G Aafl for a > 1, then / G L^(M) and by the 
dominated convergence theorem, / G C'(M). The index q characterizes the behavior near k = 0: 
a function which behaves like \k\'^ around zero is in the space Aa,q- The space Ta,q includes some 
characterization of the derivative with respect to k, which is needed in order to characterized the 
behavior near /c = 0 as shown in the next lemma. 


Lemma 5. For all f G Ta,q, the funetion 

g{k) = f{k) - 
where x is a smooth eut-off funetion with 


T<?l-i 


E 7T-«U{0) x{l»:|) 


i=0 


X([0,1]) = {1} , x([2,oo)) = {0} , 

satisfies g G Wa,p and therefore 

= {/ G ra,q: = 0, Vi < LgJ - 1 } . 

Proof. Due to the fact that the behavior at large |/c| of functions in Ta,q and in y\ia,q are the 
same we only need to prove the behavior for small |A:|. In view of the properties of the cut-off 
function, for |/c| < 1 and i < [gj — 1 , we have 


dk9{k) 




y 


dlf{k) - I I x{\k\) 

\k\ 


j=o a- 


■) 


< 


Jo 


\k\ 


q-i 


□ 
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For functions of k and y, we define the following spaces with norms reflecting the scaling 
property of the Jeffery-Hamel solution: 


Definition 6 (function spaces on fl). For a > 0 and g G M, let Ba,q be the Banach space of 

functions f G C{X x [1, oo), C) where X is defined by (6), such that f{k, y) = f{—k, y) and such 
that the norm 

\\f-Ba,q\\= sup 

Xx[l;oo) k'a,q 


is finite, where the weight is given by 


k-a,q{k,y) 


y<il + {\k\yr' 

1 1 / 1 
2/9 1 + (|A:|y)" (|A:|//)“'' 


q>o, 
q <0. 


For a > 0 and g > 0, the space for the velocity field Ua,q is Banach space of functions in Ba+i,q-i 
such that the following norm is finite, 


Bl 1«1 

-h9-i-*+ill • 

2=0 j=0 

For a > 1 and q > 1, the spaces for the pressure Va,q and for the nonlinearity TZa,q are the 
Banach spaces of functions in such that the respective norms are finite, 

[gj-l [oj 

2=0 j = 0 
[qJ-1 H-l 

\\f-,n^,q\\= E E ||9^9l/;Sa+l-,,g-l-2+,||. 

2=0 j = 0 

Remark 7. The parameter a captures the behavior of functions at infinity as a function of \k\ y, 
which is reminiscent of the scaling properties in x/y of the Jeffery-Hamel solution. By taking the 
inverse Fourier transform, the parameter a corresponds to the regularity in x in direct space. The 
index q determines the decay in y at inhnity. As we will see below, functions on the boundary 
which are in Aa,q are in the space Ba,q, when evolved in time by The spaces Ua,q, 'Pa,q 

and TZa^q include derivatives with respect to k in order to catch the behavior near k = 0 and 
derivatives with respect to y for the regularity in the //-direction. 

Remark 8. For a' > a and q' > q we have the inclusion Aq,/ g/ C Xa^q for X = A, B, T, W, U, V, 
and TZ, which will be routinely used without mention. 

Remark 9. Since the completion dehning B^^q is dehned by starting from smooth functions on 
the closed set M x [l,oo), the restriction of a function / G Ba^q to the boundary y = 1, is a 
function in Aa,q- In the same way, the restriction of / G TZa^q is in Ta,q- 
These spaces lead to the following regularity in direct space: 

Lemma 10. For a > 1 and q > 0, if f G Ba,q, we have 

fGC{n), y'+'?/GL-(n) , yi-^fGL\n), 


for all e > 0. For a and g > 0, if f G Ua,q or f G Fa,q we have 

f G (n ), 9^/ G L-(fl), e L^i^), 

for i + j < [aj, and e > 0. 
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Proof. We consider / G 13a,q. At fixed y, /{•,y) G so / is continuous in x. The continuity 

in y follows from the fact that f{k, •) G (^([IjOo)), so / G (7(0). Since 


J~[P‘a,q\ I < 


1 

yg+1 


1 , _^ 

1 + z" “ a-1 yi+^ ’ 


then G L°°. Finally, by Parseval identity 


|/(x,y)7 dx= \f{k,y)\ dk < 


\f;B. 


a,q \ 


.,2q-\-l 


l + z' 


2a 


dz< 


,2g+l 


SO y'^ ^ f G L^(0) for all e > 0. 

Finally, we consider / G Ua,q or / G 'Pa,q- Since 9*0^/ = T'[(iA:)* 5^/] and d0 G Ba-j,q+j, 
we have G Ba-j-i,q+j+i, so by applying the previous result, we obtain the claimed 

properties. □ 


3. Stokes system 

In this section we consider the following inhomogeneous Stokes system, 

Au-Vp = V - Q, V-u = 0, u\g^ = u*, (7) 

where Q is a given symmetric tensor. For simplicity we define 

R = Qi2 = Q21, ^ ~ 2 ~ ' 

The aim is to determine the compatibility conditions on the boundary data u* and on the 
inhomogeneous term R = {R, S), such that (7) admits an (a, g')-solution: 

Definition 11 ((a, g)-solutions for Stokes). A pair {u,p) G h(a,q x 'Pa-i,q+i is called an {a,q)- 
solution of the Stokes equation, if it satisfies the Fourier transform (with respect to x) of the 
Stokes equation (7). 

Lemma 12 ((a, g)-solutions are classical solutions). For a > 2 and q > 0, if {u,p) is an 
{a, q)-solution, then its inverse Fourier transform {u,p) has the regularity {u,p) G (7^(0) x (7^(11) 
and satisfies the Stokes equation (7) in the classical sense. 

Proof. In view of lemma 10, we obtain that {u,p) G (7^(11) x (7^(11), and since {u,p) satisfies 
the Fourier transform of (7), we obtain that {u,p) is a solution of (7) in the classical sense. □ 

Theorem 13 (existence of (a, g)-solutions for Stokes). We have: 

1. For all a > 2 and q > I, if Q = 0 and u* G Wa,q, there exists an {a, q)-solution. 

2. For all a > 2 and q > 1 with q ^ N, if Q & lZa,q+i and u* G Ta,q there exists an (a, q)- 
solution provided the compatibility condition Ur G y\ia,q holds, which is explicitly written in 
the following proposition. 

Proposition 14 (compatibility conditions). For q G (1,2), there is one compatibility condition 

/ CO 

i?(0,y)dy = 0, r(0) = 0. 
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and for p G (2,3), we have the additional conditions, 

/ OQ roo 

dkR{t), y) dy-2i J {y - 1 ) 5(0, y) dy = 0 , 

/ oo roo 

dkRiO, y)dy + ij {y - 1) R{0, y)dy = 0. 

The rest of this section is devoted to the proof of the existence of (a, ( 7 )-solution for Stokes 
system. 

Definition 15. We define the following operators: 

{T^w) {k, y) = ll" (1 _ ^ _ 1 ))) dz , 

{T^w) {k, y) = 2 J ±x{\k\{z- 1 )) w{k, z) dz , 

(Urw) {k, y) = {y- 1)'’ , 

where x ^ C'“(M) is a smooth cut-off function such that 

X([0,1]) = {1}, x([2,oo)) = {0}, ( 8 ) 

and their combinations: 


T+ = r< - r+, 


B+ = r+ 


y=l 


T = crr< -I- (tT^ , 
B- = T-\ , . 

\y=l 


Proposition 16 . Formally, the Fourier transform of the Stokes system is given by 
u = Ni? -|- BUj., Uj. = u* — N|y=i R ) 


p = —ikT^R -|- \kT S — Qi 2 — Uq [ 2 ik {ur -|- crvr)] , 


(9) 

( 10 ) 


where 


N = 


T+ -T-ik{z-y) -T- -T+ik{z-yy 


T~^ik {z — y) 
Proof. The vorticity is 


-T \k{z-y) 


B = 


Uq -\k\Ui -ikUl 
-ikUi UQ + \k\Ui 


uo = dxV — dyU , 


and the Stokes equation (7) implies the vorticity equation 

Acj = — dy^ R — 2 dxdyS. 


( 11 ) 


( 12 ) 


By defining 7 = a; -|- i?, the divergence-free condition, (11), and (12) can be rewritten as a first 
order differential system in y, 


dyU = dxV — x + R, 

dyV = -dxU , 


dyX = dxV - “^dxS 
dyf] = -dxX + 2<9x-R • 


By taking formally the Fourier transform in the variable x, the divergence-free condition, (11) 
and ( 12 ) can be rewritten as a dynamical system where y plays the role of time. 


dyV = Lr + q, 


u{k, 1 ) =Ub, 
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where 


r = 


h 


(° 

iA: 

-1 

0^ 


/ R \ 

V 


—ik 

0 

0 

0 

, q = 

0 

7 

? L = 

0 

0 

0 

iA; 

-2iA;5 

\v) 


1 0 

0 

—iA: 



iy 2iA:.R / 


The eigenvalues of L are given by ± |A:|, and since we are interested in solutions with zero velocity 
at infinity, we have to distinguish between stable and unstable modes, so the solution is given by 

/ y ^cxD 

z)dz — J (1 — P) z) dz + e^^rs{k) 

/ OO 

(1 - P) z) dz + ffe(fc) , 

where P is the projection onto stable modes and r* is such that the boundary condition in (7) is 
satisfied. 


P = 



(\k\ 

—iA: 

1/2 

0 ^ 


( Ub \ 

1 

\k 

|A:| 

0 

-1/2 

, h = 

Vb 

m 

0 

0 

\k\ 

—iA: 

2 \k\ {ub + avb) 


lo 

0 

\k 

\k\ j 


\ 2\k {ub + crvb) J 


where 


cj = i • sign(A:) . 


By using the Jordan decomposition for L we can explicitly calculate the exponential and we find 
that 


1 n 


1 


“ 2 


-my-z) _ 1^1 _ / ^-\k\{z-y) ^ 1^1 (y _ 


+ e-l^lto-1) [(1 - \k\ (y - 1)) Us -ik{y- 1) u,] . 

1 
2 


1 ry 1 

v = -- {y - z) R-{k, z)dz+- / {y - z) R+(k, z) dz 


+ e [(1 + |A;| (y - 1)) - iA: (y - 1) Us] , 

where R± = Rdz aS, and 


1 


Us = u + - 


^-|fc|(^-i) /I _ 


” +2 


2 
1 


{1 — \k\{z — 1)) i?+(A:, z) dz . 


'1 


{z - 1) R+{k, z) dz. 

By using the operators defined in definition 15, we can rewrite the integral equations as 


1 


u = T+ 
D = r+ 

where 


R — \k{z — y) S 


-T' 


S + \k{z - y) R +U[{1- |A:| (y - 1)) Ur -\k{y - 1) Vr] , 


iA: {z — y) R 


-T- 


ik{z-y)S +U [(1 + \k\ (y - 1)) Vr -ik{y- 1) Ur] , 


Ur = U* — 

Vr = V* - B+ 


R — ik{z — 1) S 


+ B- 


S + '^k{z — 1) R 


\k {z — 1) R 


+ B- 


\k {z — 1) S 


which shows (9). Finally, from the Fourier transform of the Stokes equation (7), we can check 
that the pressure is effectively given by (10). □ 


10 





























In order to prove the existence of an (a, g)-solution, we have to estimate the operators used in 
proposition 16: 

Lemma 17. For a > 1, r G N, and p G M, the operator Ur : Aa,q Ba+r,q-r is well-defined 
and continuous. 


Proof. It suffices to prove that 


1 


i + \kr V y 


y- 1 




1 


l + {\k\y) 


a+r * 


For \k\y < 1 the result is trivial and for \k\y > 1, we distinguish two cases: for y > 2, we have 


1 + 


1_ f y _g-|fc|(y-i) < Q-\\k\y < 


1 


1 + {\k\y) 


Oifi-r ’ 


and for 1 < y < 2, 


1 


l + \kf 


y- 1 




l + i\k\yri\k\yy 




1 


l + (|fc|y) 


a+r * 


□ 


Lemma 18. For a > 1, r G N and q >0, the operator Ur: Wa,g hla-\-r,q-r is well-defined and 
continuous. 

Proof. First of all, we have 

dyUr — Ur—\ “t“ I A: I Ur , 

so dyUr : ++i,q_i —)■ Ba+r,q-r- For (7 > 0, we have 

dk {UrW) = UrdkW + iaUr+lW , 

so by lemma 17 we obtain that dk (Urw) G Ba+r+i,q-r- 2 - The result now follows by a recursion 
on the number of derivatives. □ 

Lemma 19. For a > 0 and q > 0, the operator r<: Ba,q —)• Ba+i,q-i is well-defined and 
continuous. 

Proof. Due to the cut-off function, the integral vanishes for \k\ {y — 1) < 1, and for | A:| y > 1 -|- | A:|, 
we split the integral: 

y + 1 

(T<ya,g) {k, y) ^ e-\k\iy-z) (1 _ ^ (|/j| _ 1 ))) ^a,q{k, z) dz -h y,a,q{k, z) dz 

<^-my-i)/2 ^ r^^^fik)dz + yo.Ak,y) dz 

< Q-\k\iy-l)/2 (/j) + y) < fik, y) , 

|A:|y 

where for the last step we apply lemma 17. □ 

Lemma 20. For all a > 1 and q > d with q 1, the operators r> : Ba,q —>■ Ba+i,q-i are 
well-defined and continuous. 
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Proof. For \k\y > 1, we have 




roo PC 

{Tyfia,q) (k, y)< z) d 2 : < fj-a,q {k, y) / 

J V J V 


< 

rs^ 


1 


I; I e'a,? 

\k\y 

and for |A;| y < 1, since p ^ q, we have 


Pa,q-l{k,y) < pa+l,q-l{k,y) 


{T>Pa,q) {k,y) 


1 1 


</ Pa,q{k,z)dz<\k\’^ ^ I -d- " du< Pa+l,q-l{k,y) 

Jy JlklyUll+W^ 


□ 


Lemma 21. For all a > 0 and q > 0 with q ^ N, the operators T±: TZa,q+i Ua+i,q are 
well-defined and continuous. 

Proof. First, by using lemma 20, we have T±: 13a+i,q 13a-e2,q-i- By take the derivative with 

respect to y, we get 

dy {T<q) {k, y) = - x{\k\{y - 1))) q{k, y) - \k\ {T.^q) (fc, y) , 

dy (Tyq) {k,y) = (1 ±x(|A;| {y - 1))) q{k,y) + \k\ (Tfiq) {k,y ), 

so the time-derivative of the operators are 

dyT+ = ^+ ikP- , dyP- = -ikT+ , 

so dyT±: Ba+i,q —)• Ba+i^q. Since the integrand of Tfi vanishes at /c = 0, we have 

dk {T^w) = T^dkW -I- i ( 7 / - 1) aT<^w - iaf^ {z - l)w , 

dk (^> ^) = TydkW — i (y — 1) aTfiw + iaTfi {z — l)w , 

dk (aTfiw) = aTfidkW -|- i (y - 1) T^w - iT;^ {z -l)w , 

where a tilde over an operator denotes the same operator where x is replaced by x + which is 
also a cut-off function satisfying (8). Therefore, 

dk {T+w) = T+dkW -I- i (y - 1) T-w - \af- {z -l)w , 
dk (T-w) = T_dkW - i (y - 1) T+w -h ir+ {z-l)w, 

and by using the previously shown properties on the operators T±, we obtain that dk {T±w) G 
Ba+ 2 ,q- 2 - By recursion on the number of derivatives we obtain T±w G lJa+i,q. □ 

We can now apply these lemmas to prove the existence of (a, y)-solutions: 

Proof of theorem 13. By applying lemma 18, we have B: Wa,p —>■ hla,p, and therefore, if Q = 0 
and u* G yVa,q, we obtain that {u,p) G Ua^q x Va-i,q+i- By applying lemma 21, noting that 

z-y = {z-l)-{y-l) , 

and bounding each resulting term separately, we obtain that N: lZa,q+i hla,q, for q > 1 
with q ^ 'H. In view of remark 9, we have : lZa,q+i —)• Ta^q. In case Ur G Wa,g, since 

B : >Vq,,p l{a,p, we have {u,p) G Ua,q X Va-l,q+l. □ 

The deduction of the compatibility conditions is now straightforward: 

Proof of proposition If. Since Ur G Ta,p, we use the characterization of }Va,p in terms of elements 
of Ta,p provided in lemma 5. The first compatibility condition is ■Ur.(O) = 0, and the second 
dkUr{d) = 0. By explicit calculations, we obtain the claimed conditions. □ 
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4. Strong solutions to the Navier-Stokes equation 

The Navier-Stokes equation in the half-space can be written as the Stokes system (7) with 
Q = u <Si u, and we are going to look for solutions of the form u = u'^ + iii and perform a fixed 
point argument on ui G Ua,q- First of all, the Jeffery-Hamel solution (2) at fixed values of y and 
large values of ±s, where s = x/y, is 

ul{x,y) ~ , vUx,y) ~ —^/'(±f) > 

y s y s 

so that its Fourier transforms satisfies dyU% G Ba,i, so 

u% G lAoifi , (13) 

for arbitrary a > 1. In order to treat the nonlinearity uj we need the following proposition 
concerning the convolution: 

Proposition 22. For a > 1 and q > 1 the convolution *: Ua,i x Ua,q —?■ 'F.a,q+i is a continuous 
bilinear map. 

Proof. First, we show that the map =t=: Ba,p x Ba,q-i —>■ Ba,p+q is a continuous bilinear map, for 
p,q > 0. If / G Ba,p and g G Ba,q-i, f is in L“(M) and g is in L^(M) for fixed y G [1; oo), so (see 
for example Folland, 1999, Proposition 8.8) f * g G C(M). The dependence of the convolution on 
the power of y is trivial and it therefore suffices to prove that 

Fafi * hafi ~ ha,! j 

for L' G (0,1), where 

11 

For |A:| y < 1, we have 

[ hafi (^, y) hafi {k -i,y)di< [ ^^^0 (k - i, y) M 

>/ M. >/ M 

< - [ 
y JR 

< Fa,! {k,y) , 

and, for |A;| y > 1, we have that q — hafi and therefore, by splitting the integral at k/2, we 
find that 

/ ya,oi^,y) hafiik - i,y)di < ya,o{k/2,y) / ya,o{^,y)di 

M. M 

< ya,o{k/2,y)- [ /i« o(Al)d£ 

y Jr 

< Fa,! {k,y) . 

Now we consider / G h(a,i, and g G U^^q. If a > 0, we have (see for example Folland, 1999, 
Exercise 8.8) dk (^f * g^ = f* dkg, so by using the previous result, f * dkg G Ba+i,q-i. By taking 
the derivative with respect to y, we have dy (^f * g^ = dyf * g + f * dyg G Ba,q+i. Finally, by a 
recursion on the number of derivatives, we obtain that f * g G TZa,q+i- D 
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Now we can state the main theorem: 

Theorem 23 (existence of (a, g)-solutions for Navier-Stokes). For a > 2 and q G (1,2), there 
exists n > 0 such that for any G M and Ug G Ta,g satisfying 

Us{x) d® = 0 , 

there exists ^4 G M such that there exists {u,p) G C^(0) x (7^(0) satisfying (1) with 

u*{x) = u%{x, 1) + e~2^^ + Us{x). 

V 2vr 

Moreover, ii — -uj G Ua^q so that 


l^hl < z/, 


II — ^ ; 


and 


lim y sup \u — uj 

y^°° V^eiR 


= 0 , 


Vu G L‘^{n ), 
u/y G L^(n), 


yu G L“(S4), 
y'^Vu G L“(fl). 


(14) 


Proof. We look for solutions of the form u = + ui and perform a fixed point argument on 

■ui in the space ui G l'(a,q. In view of the previous section, the Navier-Stokes equation can be 
written as the Stokes equation (7) for ui where Q = 0ui + ui0 itj + ui<SiUi. The boundary 
condition is 

u*{k) = u%{k, 1) -f {A, 0) e“5^^ Us{k ), 
and the compatibility conditions of proposition 14 are given by 


'^r(O) — H“ 



RiO,y) dy, 


{)r-(0) = 0, 


since by hypothesis Us{k) = 0. Therefore, by defining A = — R(0, y) dy, the two compatibility 

conditions are fulfilled. In what follows, C > 0 represents a generic constant depending on q, but 
not on e. By proposition 22, we have Q G TZa^q+i and 


'R'a,q+l ^ C (II tt^, 77 q,^o II “f II Hi, II) l|Hl, Ua,q \ 


Since 


we have 


1^1 < 


R{d,y) dy 


< 


P 


R, Roi,q-\-l 


\\u* - ■&!; Ta,q\\ < \\Us; Ta,q\\ + C 
By applying theorem 13, we obtain that 


R: TZ. 


5 


II II ^ C \ \Ug , II H“ C (|| ; ^q,0 || “1“ ||"^1; ||) 11 *^1; || * 

Therefore, for e > 0 small enough, a fixed point argument shows the existence of a solution 
{u,p) G Ua,q X 'Pa-i,p+i of the Fourier transform of the Navier-Stokes equation. In the same way 
as in theorem 13, we obtain the claimed regularity and the asymptotic properties. □ 
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5. Existence of weak solutions 


In this section we define weak soiutions for our probiem, and we discuss in particuiar the 
technicaiities due to the inhomogeneous boundary conditions on an unbounded boundary. In 
order to show that our definition of weak soiutions is generai enough, we then construct such 
soiutions by Leray’s method. To study an inhomogeneous boundary probiem, it is standard (see 
for exampie Ladyzhenskaya, 1969, Chapter 5.) to dehne weak soiutions by using an extension 
map to write the energy inequaiity. 

We denote by the subspace of the homogeneous Soboiev space of order (1,2) of 

divergence-free functions on Q, and by Dq‘^{Q) the compietion with respect to the norm of 
Hh^(fl) of the vector space of smooth divergence-free functions with compact support in fl. We 
refer the reader to (Gaidi, 2011, Chapter II.6.) for the properties of these spaces. The main tool 
in studying the existence and uniqueness of weak solution is the Hardy inequality: 

Proposition 24 (Maz’yaj 2011, §2.7.1). For all u G Ilg’^(H), we have 

\\u/y\\2 < 2 ||Vw||2 . 

We define an extension map as foiiows, whose existence is proved in theorem 13 for ii* G Aa,o- 

Definition 25 (extension). Given a boundary condition u*, an extension is a map a G Zl^’^(H) 
such that a/y G L^(H), ya G L“(H) and y'^Wa G L°° and such that the trace of a on dFl is u*. 

Definition 26 (weak solution). A weak solution in the domain H with boundary condition u* 
is a vector field u = a + where a is an extension of u* and v G ZIq’^(H) which satisfies: 

/ Vri:V(/?-|- / (li • Vri) • (/? = 0 , (15) 

J ^ Jq 

for arbitrary smooth divergence-free vector-fields (f with compact support in H. 

The main result of this section is the existence of weak solutions: 


Theorem 27 (existence of weak solution). For a small enough boundary condition u* (more 
precisely such that ||?/a||oo + ||a/y ||2 is small enough), there exists a weak solution u in H. 

Before proving this theorem, we mention the fact that any weak solution vanishes at infinity 
in the following sense: 

Proposition 28. If u = a + v is a weak solution with v G ZIq’^(H), and ya G then 

lim u = 0 , 

|a;|—>00 


in the following sense 


lim 

r^oo 



u{r sin 9,1 + r cos 0) d0 = 0 . 


Proof. First of all, by Hardy inequality, we have u/y G L?‘{VL). We dehne the half-ball and 
the half-shell Sn by 


Q.n = B{{0,1 ), n) n H , 


Sn — ^2n \ fin ; 


with H((0, l),n) the open ball of radius n centered at (0,1). By using the trace theorem in Si, 
there exists C > 0 such that 


u-,L?‘{dFt\) < u-,L^{dSi) 


< C 


u-,L\Si) 


+ C 


^u-L\Si) 
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By a rescaling argument, we obtain that 


- u-,L?{dVtn) <w;L^(5„) +C Vu-,L^{Sn) , 
n 

and since y < n in we have 

- u-,L?{dVtn) <C u/y,L^{Sn) +C Vu]L^{Sn) ■ 
n 

In the limit n —?■ oo, the right hand-side converges to zero, because u/y, Vri G and since 

the integrals over Sn can be written as the difference of integrals over Q 2 n and 0,n- Finally, 


|ii(r sin 0,1-|-r cos 0)1^ d0 = -— u]L^{dVtn) , 


and the result is proved. 


As usual, to show the existence of a weak solution in an unbounded domain, we first prove, 
for arbitrary n G N, the existence of a weak solution in the domains 17^ defined in the previous 
proof. To this end, we introduce the concept of approximate weak solution in and then apply 
the Leray-Schauder theorem to prove the existence of such approximate solutions. 

Definition 29 (approximate weak solution). For n G N, an approximate weak solution is a 
vector field Un = a + v^, where G Dg’^(n) with support in which satisfies 


/ • Vu„) • (^ = 0 , 

J Q, J Q, 


for arbitrary smooth divergence-free vector-fields with support in 

Lemma 30 (existence of approximate weak solution). Provided u* is small enough, there exists 
for all n G N an approximate weak solution Un = a + Vn, with || Vr^nll < 1- 

Proof. First we note that the trilinear term can be bounded as 


Therefore the map 


{Un-'VUn)-(f < ||Vlln||2 ||W„||4 ||(^| 


HU^r. 


{Un- VUn) -ip, 

JQ 

is a continuous linear form, and by the Riesz representation theorem, there exists a map 
Fn : IFo,a(^n) ^ such that 


{Fn{Vn),<p) = - {Un-VUn)-ip. 

Jn 

The map F^ is continuous on IFQ^’^(fI„) when equipped with the L^-norm, and, since IFg^’^(fI„) 
is compactly embedded in L^(n„), is completely continuous. 

The problem of finding an approximate solution is equivalent to solving the equation 

— Fn{Vn) 
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in WQ’‘^{Vln)- From the Leray-Schauder fixed point theorem (see for example Gilbarg &: Trudinger, 
1988, Theorem 11.6.) to prove the existence of an approximate weak solution it is sufficient to 
prove that the set of all possible solutions of the equation 

Vn = \Fn{Vn) , (16) 


is uniformly bounded in A E [0,1]. 

To this end, we take the scalar product of (16) with and after integrations by parts, we get 

/ Vu„ : Vu„ = A / {un-Vvn)-a. 

Jn Jn 

Therefore by Holder inequality, we obtain 

||Vu„||2 < A(||a/y||2 + Iji’n/ylla) l|Vu„||2 ||ya||^ , 

and therefore by using Hardy inequality, 

||Vu„||2 < A(||a/y||2 + 2 ||Vu„||2) ||ya||^ . 

For A E [0,1], we finally obtain for n big enough, and a small enough, 

|a/y|l2 lly«lloo 


|VU„||2 < 


l-2||ya|loo 

which proves that is uniformly bounded. 

We are now able to take the limit n —)• oo and prove the existence of a weak solution in H: 


□ 


Proof of theorem 27. By lemma 30, there exists for any n E N an approximate weak-solution 
Vn and the sequence (u„)^gj,j is bounded in ZIq’^(H). Therefore, we can extract a subsequence, 
denoted also by (u„)^gpj, which converges weakly to v in ZIq’^(H). Now let be a test function 
with compact support in H. Then, there exists m E N such that the support of (p is in Hm- 
Therefore, we have for any n > m, 

/ Vu„:V(/?+ / {Un •'Vun) • if = 0 . 

Jn Jn 

By replacing the test function by Vn and after integration by parts, we obtain 


/ VVn ■■VVn < / {Un-'VVn)-a. 

J Q J Q 


Therefore it remains to prove that these last two equations remain valid in the limit n —>■ oo. By 
definition of the weak convergence, we have 


and 


lim / VVn'.'V(f= / Vu:Vo?, 

Jn 


Vu : Vu < lim inf / : Vu„ < 1, 

, n^OO 
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Since has support in Om, 


/ {Un • VUn - U • Vu) • if 

— 

/ {{Un-U) -VUn + U 

■V {Un- U) ) • (f 

Jq 


J ^rn 




< 

/ {{Vn-V) -VUn) • if 

+ 

1 {u-V(f) 







< (||Vw„||2 llv^lL + 2 II Vunlla l|yV(^||^) ||u„ - V] L2(0^ 


and therefore since is compactly embedded in L 2 {^m) this proves that 


lim / (n„ • Vun) ■ip= / (w • Vu) • >p>. 


□ 


6. Uniqueness 


In this section, we prove a weak-strong uniqueness theorem by exploiting the properties of 
(a, ( 7 )-solutions. Namely we prove that any weak solution satisfying the decay properties (14) of 
an (a, g)-solution coincides with any weak-solutions for the same boundary data. The ideas of 
the proof that are not specific to the presence of an extension can be found in Hillairet &: Wittwer 
(2012) and we refer the reader to this article for some technical details which are omitted here. 

Theorem 31 (weak-strong uniqueness). Let he u a weak solution that satisfies 

ueDl’\n), yu€L°°{n), y‘^Vu€L°°{n), (17) 

and such that ||yu||^ is small enough. Then any weak solution u with boundary value u* = 111^^ 
that satisfies the energy inequality 


coincides with u. 


Vri : Vu < / (« • Vu)a, 




(18) 


Remark 32. The (a, g')-solutions found in theorem 23 satisfy the requirement (17) on u. 

The remaining part of this section is devoted to the proof of this theorem. To begin with, we 
prove that integration by parts with respect to the solution u is permitted: 

Lemma 33 (integration by parts). For any u that satisfies (17), we have 


/ {w-'Vu)-u+ / {w-'Vu)-u = 0. 

Jq Jq. 


for all u,w G D]p‘^{Tt) with ujy £ and w/y £ L^(fl). We note in particular, that if u and 

w are weak solutions, the hypothesis are satisfied. 

Proof. By using Holder inequality, we have the bounds 


and 


/ {w •’Vu) • u 

Jq 


< W^/vh 


y'^'Vu 


\u/y\\ 


2 ) 



{w • Vri) • u 


< \\w/y\\2 l|Vu||2 I|yw||^ 
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Then, the result follows by an integration by parts, where u is approximated by compactly 
supported functions (see for example Hillairet &: Wittwer, 2012, Proposition 21). 

Finally, if u = a + i; is a weak solution, we have by hypothesis a/y £ Lp‘{yi) and by Hardy 
inequality vjy £ since v G □ 

Next we prove some results on the extension of allowed test functions in the definition of weak 
solutions: 

Lemma 34. If u is a weak solution, then 


/ Vw : Vt; + / (it • Vit) • a = 0 , 

Jq Jfl 


for any v G sueh that yv G 

Proof. We have 

/ Vri:V(^+ / {u-'Vu)-ip < / Vu : V(^ 

Jq Jq Jq, 


+ 


(a • Vu) • (f 


+ 


{v • Vw) • (p 


< IIVu||2 + llya|loo Il^^^ll2 \\p/y\\2 

+ \\v/y\\2 llVwlla ||2/(p|loo 

< l|Vw||2 (II V(p||2 + 2 Hyalin || V¥j|| 2 + 2 ||Vi>||2 ||y<p||oo) , 

so that the expression under consideration defines a linear form in p. Since we can approximate 
V by (hn)„gpj G (^“^.(H), such that || Vh — Vhn ||2 + \\yv — yVnWao —^ 0 as n —)• oo, this proves 
the lemma. □ 

Lemma 35. If u is a weak solution sueh that yu G L°°(H), then 


/ Vit : Vr; + / {u-'Vu)-v = 0, 

J Jq. 


for any v G 
Proof. We have 


/ Vti:V(^+ / {u-'Vu)-(p < / Vu : 

Jfl Jo. Jo. 


+ 


{u • Vri) • If 


< IIVw||2 IIV(^||2 + ||yw||oo IIVw||2 \\f/y\\2 

< l|Vw||2(l + 2 ||yw||^) ||V(^||2 , 

and since the form is linear in (f, the lemma is proved. 

We now prove that the weak solution u satisfies an energy equality: 

Lemma 36. Any weak solution u whieh satisfies (17) verifies the energy equality 


□ 


Proof. By lemma 35, we have 


/ Vu : Vn = / {u-'Vv)a. 

Jn Jn 

/ Vw : Vh + / {u-'Vu)-v = 0, 
Jn Jn 


(19) 
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and by lemma 33, 


/ {u-'Vv)-v = 0, / {u-'Va)-v+ / {u-'Vv)-a = 0, 

Jo, Jo, Jfl 

so we obtain the energy equality 

/ Vii : Vh = — / {u • Va) • v = / (ii • Vh) • a . 

Jvt J Q. Jo. 

We now have the necessary tools in order to prove the main theorem of this section: 


□ 


Proof of theorem 31. Let u and u be two weak solutions with the same boundary conditions, so 
d = u — u £ see for example Galdi (2011, Theorem II.7.7). Then, by using the scalar 

product on we have 

||Vd||2= [ (Vu-Vu) : {Vv -Vv) 

Jn 

= / Vw : Vt; + / Vfi : Vh — / Vit : Vh — / Vw:Vt;. 

Jn Jn Jn Jn 

By using lemmas 34 and 35, the energy equality (19) and the energy inequality (18), we have 
l|Vd ||2 < / (u-Vtj)‘a + / (ii‘Vr;)‘a+ / (ri‘Vu)-t>+ / {u-^u)’V. 

J fi Jq Jq, J fi 

Since ujy G Lp‘{Vl) and ujy £ L^(II) by using Hardy inequality and lemma 33 we have 

/ (ri-Va)-a = 0, / ('a-Va)-a = 0, 

v/o J u 

which allows us to rewrite the bound as 

\\Vd\\l< / {u. Vu) • a+ {u- Vu) • a+ {u • Vu) • V + / (w • Vw) • V 
Jn Jn Jn Jn 

< / •Vti)'it+ / {u-Vu)-u. 

Jq Jq 

By using lemma 33, we integrate the second term by parts, 

||Vd||2 < / {u-Vu)-u— / {u‘Vu)‘U = / {d‘Vu)‘U. 

Jn Jn Jn 

Again by lemma 33, we have 

/ (d • Vti) • w = 0 , 

Jn 

so by Hardy inequality. 


l|Vd||2< 


(d • Vd) • u 


<||d/y||2||Vd||2||yd|L<2||yd||^||Vd||^ 


Therefore, if ||?/w||^ is small enough, we obtain that d = 0, i.e. u = u. 


□ 
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7. Numerical simulations 


In order to simulate this problem numerically, we truncate the domain H to a ball of radius 
R = 10^, = Oni?((0,1), -R). On the bottom boundary we take an antisymmetric perturbation 

of a symmetric Jeffery-Hamel, 

+ ^sin(20)e^, (20) 

and on the artificial boundary F, which is the upper half circle of radius R, we take 

ll|p = . 

In figure 3, we represent the velocity field u multiplied by r in order to see the behavior at large 
distances. For 1 ^ = 0, this corresponds to the Jeffery-Hamel solutions which are scale-invariant. 
For negative fluxes $ < 0, small perturbations have almost no effect on the behavior at large 
distances, so the asymptotic term is probably given by the Jeffery-Hamel solution Conversely 
for > 0, even a small perturbation drastically change the behavior of the solution at large 
distances by somehow rotating the region where the magnitude of the velocity is large. In the 
case, the asymptotic behavior is very likely not given by the Jeffery-Hamel solution This 
conclusion can also be seen in hgure 4, where we plot the velocity held u in polar coordinates 
multiplied by r on the half-circle {10^ (sin 0, cos 0) ,6 E [^; |]} in terms of v. 


21 


$ =-2,1/= 0.0 ^> = -2,1/ = 0.1 $ = -2, z/ = 0.2 



^ = +1, u = 0.0 


$ = + 1 , 1 / = 0.1 


$ = +1, z/ = 0.2 


10^ 


0 




0 io3 


$ = + 2 , 1 / = 0.0 


$ = + 2 , 1 / = 0.1 


$ = +2, z/ = 0.2 


10^ 


0 





0 io3 


Figure 3: Numerical results for the velocity u multiplied by r, for the boundary condition (20) 
in a domain of size R = 10^, for various values of the flux <1> and of the perturbation v. 
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Figure 4: Profiles of rUr and rug on the half-circle of radius 10^ in term of 0 G [—7r/2,7r/2] and 

[ 0 , 1 ], 
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A. Jefrery-Hamel solutions with small flux 

A Jeffery-Hamel solution in the upper half-plane, 

|(r sin0,r cos 0) , r > 0 and 6 G (^, |)|, 

is a radial solution of the Navier-Stokes equations with zero velocity on the boundary and whose 
velocity norm is f{6)/r. More explicitly / has to satisfy the boundary value problem (21). Here 
we prove an existence theorem: 

Theorem 37. For every small enough value of the flux f, the Jeffery-Hamel equation 

/" + /'+4/= 2C, /(^)=0, (21) 

admits a symmetrie solution, 

^ cos‘^{0) + 0((/>^/^) , 

and in addition if f <t), two quasiantisymmetric solutions, 

f}\e) = ±^l^sm{2e) + 0{f). 

Proof. First of all the solution of the linear equation 

r+4/ = 5, 


is given by 
where 


f{9) = A sin(20) -|- B cos(20) — L[g] (9) 


1 


1 


/ g(s) sin(2s) ds — - sin(20) / g{s) cos{2s) ds . 


'-7r/2 


<--kI2 


Therefore the Jeffery-Hamel equation and the boundary condition (21), can be rewritten as 

^+-kI2 


f{9) = Asm{29) + Ccos‘^{9) + L[f]{9 ), 
By defining 


/ +7r/^ 

f{9) sm{29) d9 = 0. (22) 

-■it/2 


fo{9) = Asin(2(9) , fi{9) = Ceos‘S(9) - ^cos'^((9) , 


f — fo + fi + f, 


the flux condition 


|>-|-7r/2 

/-7r/2 


f{9) d9 


directly gives the definition of C in term of the flux, 

2 f A ^ 2 /■+^/2 


C — -1— - 

TT 4 TT 


-7r/2 


m d9, 
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and the two equations (22) can be rewritten by substitution as: 


/ = l[(2/o + /i + /) (/i + /) 

A[(t>+ {a - 2fo{e) - 2/i(0) - m) m d9 . 


(23) 

(24) 


In order to find a fixed point of theses equations, we first solve the left-hand-side of the second 
equation 

^0 ° 

for ^0- This equation admits the solution Aq = 0 and in addition if (/> < 0 the two solutions 



Given one of these three solutions, we define 


A = Ao +A, 


and (24) becomes: 


AScj) + SttAq 


r+7rl2 

l-n/2 


[a - 2fo{9) - 2/i(0) - m) m d9-^ (3^0 + a) A' 


(25) 


It is easily verified that the maps defined by (25) and (23) map the ball 


= eMxC7°([-f,f]): U| <50(/>and |/| < 


into itself, provided 4> is small enough. Moreover, since the maps (25) and (23) are multilinear 
affine maps of A and /, they are contractions from B^j, into itself, for (j) small enough. The 
first order terms which we explicitly computed above prove the claimed leading terms of the 
symmetric and quasiantisymmetric solutions. □ 
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